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Abstract 



Using the ARGUS detector at the e~ storage ring DORIS II, we have measured the Michel parameters p, 
^, and for — * I'^vu decays in r-pair events produced at center of mass energies in the region of the T 
resonances. Using — > p^v as spin analyzing tags, we find pe = 0.68 ±0.04 ±0.08, = 1.12 ± 0.20 ± 0.09, 
= 0.57 ± 0.14 ± 0.07, = 0.69 ± 0.06 ± 0.08, = 1-25 ± 0.27 ± 0.14 and ^(5^, = 0.72 ± 0.18 ± 0.10. In 
addition, we report the combined ARGUS results on p, ^, and ^5 using this work und previous measurements. 



1 Introduction 

A general ansatz to describe purely leptonic r decays is the well-known formalism by Michel [|l|. 
This formalism includes extensions of the Standard Model such as scalar and tensor coupling to 
left- or righthandcd leptons. The derived differential width in the r rest frame neglecting radiative 
corrections and terms proportional to /m^ is 
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where x = 2Ei/mT is the scaled energy of the charged lepton, P-r the t polarization and 1? the 
angle between the polarization vector and the charged lepton momentum. Equation is valid for 
vanishing neutrino masses. The Standard Model with a pure V ~ A interaction predicts p = 3/4, 

= 0, ^ = 1, and = 3/4. 
Up to now, measurements of the Michel parameters ^ and £,S have been made by ARGUS, ALEPH, 
L3, SLD and CLEO |§, ||, In this paper, we present measurements of p, and ^(5 with a new 
tagging method [|| separately for electrons and muons. The results are obtained using tau pairs 
produced in e+e" t^t^ {p^i^r){Pi^i'^T)- The decay p^i^r is used as a tag, i.e. for 

the event selection and as r spin analyzer. They are statistically independent of previous ARGUS 
measurements using lepton tags |Q and oi-meson tags in its decay into three charged pions 



2 Method of the Measurement 

As can be seen from Eq. (|l|), the measurement of f and requires the knowledge of the t po- 
larization. At DORIS energies, the average polarization of r leptons is zero, hence no information 
on ^ and ^(5 can be extracted from the analysis of single r-decays. Only p and 77 are accessible 
from the analysis of the momentum spectrum. However, r-leptons are produced in pairs resulting 
in a spin-spin-correlation. At ^/s = 10 GeV, r pairs are produced via a virtual photon and both 
T spins are predominantly parallel due to the vector coupling, with only a small probability of 
(l — f3r)/2 ~ rn^/s ^ 3% for being antiparallel. In the analysis described here we use the decay 
p^Vr in order to obtain information on the r-spin orientation. The basic method is similar 
to the one used in the ARGUS analysis with a^Vr as spin analyzer and is described in 

more detail there. 

The performance of various analysis methods has been studied using a Monte Carlo simulation. 
Tau pair events including radiative photons were produced by the KoralB/Tauola generator [Q. 
These events passed the ARGUS detector simulation Simarg [D and the reconstuction program 
Arg13 I^, and were finally submitted to the same selection procedure as the data, which is de- 
scribed in section ^ below. 

The parallel orientation of the two r spins is used to relate the spin of the tau decaying into l^vD 
to the spin of the tau decaying into p^i^r- The spin in the latter decay is related to the p^ direction 
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because of parity violation with a neutrino helicity h,y^ = — 1 as supported by measurements of the 
sign in —> a^Vr decays |^ and the modulus in r+r^ {p^Dr){p~i'T) events The p- meson 
can be produced in two helicity states, h = and h ~ +1. The state with h — dominates by a 
factor of (1 + — = m'^/rn^ w 5.3, where (3* is the velocity of the r in the p rest frame. 

Therefore, its r-spin p-momcntum correlation also dominates, resulting in correlations between the 
momenta pi of the lepton and Pp of the p- meson of the same tau pair event. In a Monte Carlo sample 
with the Standard Model value f = 1 on the lepton side and h^^ = — 1 on the p side, we obtain a 
correlation coefficent r = Cov{pi,pp)/{a{pi) ■ <j{pp)) — (— 6.6 ± 0.4)% whereas parity conservation 
(h^^ • ^ = 0) leads to r = (—2.2 ± 0.4)% which is non-zero because of radiation in the initial e+e~ 
state. A measurement of r is, therefore, a measurement of h^^ ■ ^. 

Additional information on ^ and ^(5 is obtained by (i) using the momentum-momentum correlations 
of the vectors instead of just their absolute values, and (ii) including information on the p-meson 
helicity, since the suppressed h — +1 state results in the opposite correlation compared to the h = 
state. A discrimination between both states is possible via the decay angle for p^ — > tt^tt", since 
the distributions for both helicities are different. 

Maximal information on p, ^, and ^(5 is obtained by a single entry maximum likelihood fit using 
the selected {l^i'iy){7T^Tr^i') events with all their 10 measured quantities, the e+e^ cms energy -y/s, 
and the momentum vectors of the three observed particles, p{P), p(7i'^), and p{tt^)- The matrix 
element for the differential cross section e+e~ — > r+r" {l^i/D){'K^Tr'^i') can be written as 

\M\' = HA[L,+pL2+ VL3] + K^KC^^l^L'^f, + ^SL'^^l (2) 

where the first part is the spin-averaged contribution and the second part describes the t^t^- 
spin correlations. H and H' describe the decays p^u and p^ — > tt^tt^. Li, L2, L^, L[, 

and L'2 describe — > l^i^D decays in the most general way with zero neutrino masses. A is 
the spin-averaged matrix element for r-pair production and C the spin-spin correlation matrix 
{a,P = 1 ... 4). includes radiative corrections which are treated in a factorizable way. 

The likelihood function to be maximized contains three variables 9 = (p, ^, £,S) and N ■ 10 measured 
quantities = {si,pi{P),pi{TT'^),pi{n^j), where i = 1...N, and is the number of selected 
events, 

N N 

m = l[u{e\m,) = l[f{m,\e). (3) 

i=l i=l 

Since data were taken at different cms energies, the probability density / is taken as a function 
of scaled momenta pi/y/s. It is the projection of the full density f{m,u\9) where u denotes the 
unobserved kinematical quantities which include neutrino momenta as well as photon momenta of 
initial and final state radiation. Furthermore, the observables rrii are obtained as a convolution of 
the true values ti of the same kinematical quantities and the detector resolution function TZ. The 
final density is given by 

f{m\9)^r]{m)-M{9) J \Mit,u\9)\^ ■ rit,u) ■ n{m\t) dt du, (4) 

where V is the phase space distribution, 77(771) is the detector acceptance, and M{9) is the normal- 
ization factor given by 

Afi9) ■ J r]{m) ■ \Mit,u\9)\^ ■ V{t,u) ■ n{m\t) dtdudm = 1. (5) 

Radiative corrections including initial state radiation and internal bremsstrahlung in r ^ evv are 
included in M.(t,u\9). External bremsstrahlung of the electron is included in the resolution function 
'R-{m\t). For the application as a likelihood function, f{m\9) is rewritten in a more appropriate 
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form. The factorizing radiative corrections allow to write u|6')p • 'P{t,u) = j=o 

with 9o = 1. The likelihood of a single event i is then 

3 

J2 dj J Fj{t, u) TZ{m,\t) dtdu 

He\m,) = , (6) 

I 'ni^) Pj (^i Ti-{'m\t) dt du dm 

where the acceptance factor has been omitted in the numerator since it depends only on the ob- 
served quantities m. The integration over the neutrino momenta, which are part of the unobserved 
quantities u, is performed to a large extent analytically, leaving only an uncertainty in the rela- 
tive angular orientation of the two tau leptons. The integration over t and the remaining seven 
unobserved variables u is done by a Monte Carlo method in which for each event 450 kinematical 
configurations are tested for compatibility with the observed set of observables rn^. The seven u 
variables in the Monte Carlo integration are the momentum vector of the initial state photon, the 
r orientation angle, and the momentum vector of the photon in t — > evv^2- A configuration is 
"successful" if the set of (mi, i, u) is compatible with the kinematics of e+e^ "fi{-K^-K^ v){r^ uD^2) 
events, i.e. if the undetermined neutrino momenta can have physical values. 

These 450 tries lead to Uhit successful matches for each event, and the distribution of Uhit is shown 
in figure ^ The high value of {nhu) shows the effectiveness of the numerical integration chosen, 
and Uhit is large enough for most events to compute Li{6\mi) with sufficient precision. Events with 
Uhit < 25 are rejected in the data selection. The good agreement between N{nhit) for selected data 
and for accepted Monte Carlo events supports the validity of the assumptions used. 

The integral in the denominator of Eq. (^, for each set of parameters 6*, is also obtained with the 
help of a Monte Carlo integration which determines the four separate integrals used in the sum. 

The validity of the method has been extensively tested by using Monte Carlo events as pseudo data 
for determining the parameter p, ^, and ^(5. Even using pseudo data samples 10 times larger than 
the real data samples, we did not detect any significant biases in the determined parameter values. 

3 Data Selection 

The measurements presented here were performed with the ARGUS detector at the e+e^-storage 
ring DORIS II. The detector and its trigger requirements are described elsewhere §. The data 
sample used was collected between 1983 and 1989 in the center of mass energy region of the T 
resonances. The integrated luminosity used in this analysis is 289 pb^^ corresponding to 265 000 t 
pairs produced. 

Event selection starts with requiring two oppositely charged tracks forming a vertex in the inter- 
action region and an opening angle larger than 80° at the vertex point. Each track must have a 
transverse momentum above 150 MeV/c and point into the barrel region. We demand one particle 
to be identified as a Icpton. The other track has to be accompagnied by one or two neutral clusters 
in its hemisphere with a minimum angle of 10° to the track. A neutral cluster is defined as energy 
deposition above 100 MeV in the calorimeter. No additional cluster is allowed, neither in the lepton 
hemisphere nor in the 10° cone around the tt* candidate. In the single cluster case, accounting for 
photons from the tt*' decay merging into one cluster, we demand a minimum energy deposition of 
1 GeV. In the two cluster case, both photons are combined to form a 7r°. The two photon system 
must have an invariant mass within 100 MeV of the nominal tt*^ mass and < 9 for the tt^ mass 
hypothesis. The lepton identification was done by a likelihood method using specific energy loss, 
time of flight, amount and lateral spread of energy deposition in the calorimeter and hits in the 
muon chambers. More details on this likelihood method arc found in refs. |ll|, |l^. To ensure a good 
lepton identification, a momentum above 0.8 GeV/c for electrons and above 1.1 GeV/c for muons 
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bars represent the data. The hatched histogram shows the expectation from the Ko- 
ralB/Tauola Monte Carlo for the decay {e^iyi'){n^n'^u) and other contributing r 
channels. The solid line marks Uhu = 25, the minimum number for selected events. 
Events with a hard initial state photon lead mostly to rihu < 25. 



as well as an electron or muon likelihood higher than 0.8 is required. The charged pion candidate 
is required to have a pion likelihood above 0.6 using specific energy loss in the drift chamber and 
time-of-flight 

To reject background, we apply the following three additional event cuts. The invariant mass m^±^o 
has to be smaller than the r lepton mass. This cut suppresses gg-events. The angle between the 
lepton momentum Pi and the p-meson momentum is not larger than 176°, which rejects radiative 
Bhabha and fi pair events. The relation fc ■ [m^nissl ^ ~ /o)^ < Vtl \fs with /c = 5 and jo = 0.45 
for the single cluster case, and fc — and Jo = 0.4 for the two cluster case has to be fulfilled. In 
here, the missing mass rrimiss uses only charged tracks and the vector sum of transverse momenta 
pt uses charged and neutral particles. This cut removes 77 and radiative Bhabha events. 

After these cuts, a data sample of 3176 candidates for {e^vv){-K^'K'^v) and 2099 candidates for 
{ijl^vv){'k^-k'^i') remains. As seen in figure ^, a good agreement for lepton momenta is found 
comparing data and Monte Carlo events. However, the data selection criteria do not allow to 
remove all background sources. The selected sample still contains 10% {■k^'k^'k^v){1^vv) decays, 
3.5% {K^'it'^v){1^vv) decays, ~ 1% misidentified leptons and around 2% background from non-r 
channels, as determined by Monte Carlo. Figure || shows the tt^tt^ invariant mass distribution with 
good agreement between data and Monte Carlo events. Tau decay fractions for the background 
estimation are taken from ref . fl^ . 

4 Data Analysis 

A first information on ^ is obtained from the pi-pp correlation coefficient. The selected data sample 
has r — —(6.1 ± 1.2)% leading to ^ = 0.89 ± 0.27 without background corrections. The finally 
used likelihood function in Eq. can conveniently include background from other semihadronic t 
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Figure 2: 

Spectra of scaled momentum p;* — pi- 10 GeV/y^ for leptons. Data are shown as dots 
with error bars and accepted Monte Carlo events, generated with p = 3/4, as grey 
shaded area. 



decays in the data sample and from lepton misidentification. The modified likelihood function is 

L* = (f - Aqi - Xk*) (^(1 - A/) • Li + Xf ■ L{°''"'^ + Aqi • Li^(^±7r07r0) + Aif. • Li^i^K±ivO), (7) 

where Aai is the fraction of the (7r^7r''7r'')-tags, Xk' is the fraction of the (i^^7r°)-tags and Xf consid- 
ers pions misidentified as leptons. Table ^ gives the used values for Xf, A^j, and Xk* as determined 
by Monte Carlo simulation. The additional likelihood functions Li k for k = {■n^'K^T:^v){r^vi') 
and {K^'K'^v){r^vD) are determined by the same technique as Li for (t:^'k'^v){1^i'D). j^i^"-^'^ jg 
the contribution to the likelihood function for misidentified leptons. Non-r background has not 
been considered in the likelihood function, because its contribution is small and is included in the 
systematic error studies, see Table ^. 









A/ 

Aai 

Xk- 


(0.8 ±0.2)% 

(10.1±^:«)% 
(3.7 ± 1.0)% 


(1.6 ±0.2)% 
(10.5 ± ^-9)% 
(3.4 ± 1.0)% 



Monte Carlo determined background levels as used in the likelihood function Eq. (4). 
The asymmetric errors are due to the systematic deviation of ARGUS branching frac- 
tions from the world average used to calculate these factors. 

The maximum likelihood fits lead to results and statistical errors for the Michel parameters as 
shown in Table ^ below. Figures and ^ illustrate the quality of the fit and the validity of the 
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spectrum includes background from other tau decays. 



assumptions made for the determination of the kinematical quantities used to calculate the matrix 
elements. 

Figure |l| shows the distribution of successful tries to find a matching initial state to the observed 
final state comparing data and Monte Carlo events. Figure ^ shows the agreement in — 21nL* 
for the fit to the data compared to the fit to Monte Carlo events generated with Standard Model 
parameters described with Eq. (0). Deviations are found at large values of —2 • InL* if we neglect 
other T sources and misidentified leptons. 

The systematic errors are summarized in Table || including limited Monte Carlo statistics, trigger 
efficiencies, branching fractions, particle identification efficiencies, lepton misidentification, the in- 
fluence of background from non-r physics and the center of mass energy dependence of the matrix 
element description. The total error is obtained as a quadratic sum of all errors in Table |^. 

5 Results and Conclusions 

The Michel parameters as found by the maximum likelihood analysis described above, and including 
systematic errors, are given in Table ^j. 

The results are in agreement with previous measurements and the Standard Model predictions. In 
addition, the errors are comparable to the errors of previous measurements. The results with p 
constrained to the world average of p = 0.742 ± 0.027 ||l^ are very close to the unconstrained ones; 
they are given in Table for completeness. 
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Figure 4: 

The contribution of single events to — 21nLj for events with rihit > 25. The points 
with error bars represent the data. The open triangles show the expectation from the 
KoralB/Tauola Monte Carlo for the decay {e^vv){Ti^-K°v) and other contributing 
r channels. The shaded area are Monte Carlo events without consideration of any 
background from other r decay channels or lepton misidentification. 

All quoted values for ^ and ^(5 assume a r neutrino helicity h^^ = —1. Using the experimental 
average = -1.011 ± 0.027 11 would not change the results. 

Assuming lepton universality and combining the results of this analysis in Table ^ with previous 
ARGUS measurements jl^ of p, ^, and ^(5, we obtain 

p = 0.731 ±0.031 
^ = 1.03 ±0.11 
S,5 = 0.63 ±0.09 

where statistical and systematical errors are combined in quadrature. Correlations in the systematic 
errors of different measurements are taken into account. The covariance matrix is given by 

p 9.54- 10--* 

C -2.30-10-5 1.25-10-2 

-7.88-10-5 1.03-10-3 7.84-10-3 
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systematic error source 


(e± 


vv) (tt^ 




{^^^ 


vv) (tt^ 






Mp) 


vs/ 










Monte Carlo statistics 


ibO.Of 


±0.05 


±0.03 


±0.02 


±0.07 


±0.04 


Trigger efficiency 


±0.01 


±0.01 


±0.01 


« 


±0.07 


±0.01 


Branching fractions 


+0.020 
-0.003 


+0.041 
-0.007 


+0.024 
-0.004 


+0.026 
-0.004 


+0.044 
-0.007 


+0.057 
-0.009 


energy dependence 


±0.03 


±0.02 


±0.02 


±0.03 


±0.03 


±0.02 


particle identification 


+0.05 
-0.06 


±0.03 


±0.03 


±0.04 


±0.02 


±0.02 


lepton misidententification 


0.005 


±0.01 


±0.004 


0.005 


±0.02 


±0.005 


background (non r physics) 


±0.01 


±0.03 


±0.03 


±0.01 


±0.07 


±0.06 


tot al systematic (>rror 


±0.08 


±0.09 


±0.07 


±0.08 


±0.14 


±0.10 



Table 2: Contributions to the systematic error 





(7rT7rOi/)(e±z/i7) 




combined fit 


p 


0.68 ±0.04 ±0.07 


0.69 ±0.06 ±0.06 


0.68 ±0.04 ±0.07 




1.11 ±0.20 ±0.08 


1.26 ±0.27 ±0.14 


1.18 ±0.16 ±0.10 




0.56 ±0.14 ±0.06 


0.73 ±0.18 ±0.10 


0.63 ±0.11 ±0.08 



Table 3: Fitted Michel parameters 
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